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Abstract 



We present a model-independent calculation of hadron matrix elements for 

all dimension-six operators associated with baryon number violating processes 

using lattice QCD. The calculation is performed with the Wilson quark action 

in the quenched approximation at (3 = 6/g 2 = 6.0 on a 28 2 x 48 x 80 lattice. 

Our results cover all the matrix elements required to estimate the partial 

lifetimes of (proton, neutron)^(-7r, K, rj) +(i>, e + , fi + ) decay modes. We point 

out the necessity of disentangling two form factors that contribute to the 

matrix element; previous calculations did not make the separation, which 

led to an underestimate of the physical matrix elements. With a correct 

separation, we find that the matrix elements have values 3 — 5 times larger 

than the smallest estimates employed in phenomenological analyses of the 

nucleon decays, which could give strong constraints on several GUT models. 

We also find that the values of the matrix elements are comparable with the 

tree-level predictions of chiral lagrangian. 
11.15.Ha,12.38.Gc,13.30.-a 
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Typeset using REVTf^X 



I. INTRODUCTION 



Nucleon decay is one of the most exciting predictions of grand unified theories (GUTs) 
regardless of the existence of supersymmetry (SUSY). Although none of the decay modes 
have been detected up to now, experimental efforts over the years have pushed the lower limit 
on the partial lifetimes of the nucleon. Moreover, an improvement by an order of magnitude 
is expected from the Super-Kamiokande experiment, which can give a strong constraint on 
(SUSY-) GUTs. On the other hand, theoretical predictions of the nucleon partial lifetimes 
suffer from various uncertainties. One of the main sources of uncertainties is found in the 
evaluation of the hadron matrix elements for the nucleon decays (PS\(D$\N), where PS 
and N denote the pseudoscalar meson and the nucleon, and 0$ is the baryon number 
violating operator that appears in the low-energy effective lagrangian of (SUSY-) GUTs. 
The matrix elements have been estimated by employing various QCD models. Their results, 
however, scatter over the range whose minimum and maximum values differ by a factor of 
ten Therefore, a precise determination of the nucleon decay matrix elements from the 
first principles using lattice QCD is of extreme importance. 

In lattice QCD the pioneering studies for the nucleon decay matrix elements 0|| at- 
tempted to estimate the matrix element (tt \O^\p), which is relevant to the dominant decay 
mode p — > 7T° + e + in the minimal SU(5) GUT, from the matrix element (0\O^\p) with the 
aid of chiral perturbation theory. This was followed by a direct measurement of (n \O^\p) 
with the use of the three-point functions [|J, which showed an unexpectedly large discrep- 
ancy between these two methods: the direct method yielded a value of the matrix element 
two or three times smaller than the value obtained by the indirect method. Recently we 
have revisited this old problem and confirmed || the peculiar feature when one follows the 
methods employed in the earlier work |4[|. 

In this paper we report results of our effort to advance the lattice QCD calculation of the 
nucleon decay matrix elements in several directions. We point out that there are two form 
factors that contribute to the matrix element (ir \O^\p) for general lepton momentum. While 
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only one of the form factors is relevant for the physical amplitudes as the other form factor 
contribution is annulled by the negligibly small lepton mass, the two contributions have to 
be disentangled in the lattice QCD calculation. This explains the discrepancy between the 
direct and indirect estimations of the proton decay matrix element found in the previous 
studies HJ^] where the separation was not made. 

Another important feature of our calculation is model independence. All dimension-six 
operators associated with baryon number violating processes are classified into four types un- 
der the requirement of SU(3)xSU(2)xU(l) invariance at low-energy scales f|[|]]. If one spec- 
ifies the decay processes of interest, namely the processes among (proton,neutron)^(7r, K, rf) 
+(u, e + , we can list a complete set of independent matrix elements in QCD, and we 
calculate all the matrix elements. 

Other advances, which are more technical but essential for precise calculation, are the 
following two points: (i) flavor SU(3) breaking effect in the process with the K meson in the 
final state is correctly taken into account by setting the strange quark mass non-degenerate 
with the up and down quark mass, (ii) two spatial momenta are injected to investigate the 
q 2 dependence of the matrix elements, where q is the four-momentum transfer between the 
nucleon and the pseudoscalar meson. 

This paper is organized as follows. In Sec. |T|we formulate our calculational method of the 
nucleon decay matrix elements. The complete set of the independent matrix elements is also 



presented. In Sec. |T| we briefly review the chiral lagrangian for the baryon number violating 
interactions and enumerate its tree-level predictions. Section [IV| contains the simulation 
parameters and technical details. Results for the matrix element (0\O^\p) are given in 
Sec. |V| In Sec. [VI] we present the results for the nucleon decay matrix elements obtained by 
the direct method and compare them with the tree-level predictions of the chiral lagrangian. 
We also discuss the soft pion limit of the matrix elements. Our conclusions are summarized 
in Sec. |VTJ. 
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II. FORMULATION OF THE METHOD 
A. Independent matrix elements for nucleon decays 

One of the most important features in the study of the baryon number violating processes 
is that the low energy effective theory is described in terms of SU(3)xSU(2)xU(l) gauge 
symmetry based on the strong and the electroweak interactions, which enables us to make 
a model independent analysis. Our interest is focused on the dimension-six operators which 
are the lowest dimensional operators in the low energy effective Hamiltonian: operators 
associated with the baryon number violating processes must contain at least three quark 
fields to form SU(3) color singlet, and then an additional lepton field is required to construct 
a Lorentz scalar. Higher-dimensional operators are suppressed by inverse powers of heavy 
particle mass that is characterized by the theory beyond the standard model. 

All dimension-six operators are classified into the four types under the requirement of 
SU(3)xSU(2)xU(l) invariance 00): 



o (1) - 

w abcd — 


( D iaR U jbR)(Q C a kcL L PdL)eijke a p, 


(1) 


o {2) - 

w abcd — 


{Q C aiahQ 


PjbL) {U% cR L dR )e ijk e al3 , 


(2) 


e> (3) - 

^abcd — 


(QaiaLQ 


PjbL) (Q^kcL^SdL)^ijk e a8^P'r, 


(3) 


^abcd — 


( D iaRUjbR) (U^ cR L dR )e ijk , 


(4) 



where %p c = ip T C with C = 7472 the charge conjugation matrix; i, j and k are SU(3) color 
indices; a, /3, 7 and 5 are SU(2) indices; a, b, c and d are generation indices; Ll and Ql are 
generic lepton and quark SU(2) doublets with the left-handed projection P L — (1 — 7s)/2; 
Lr, Ur, and Dr are generic charged lepton and quark SU(2) singlets with the right-handed 
projection Pr = (l+7s)/2. Fierz transformations are used to eliminate all vector and tensor 
Dirac structures in eqs. (0) — 

The operators relevant to non-strange final states are [|J 

°d ] = (diRUjR){u c kL e dL - d c kL v dL )e ijk , (5) 
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Of = {di L Uj L ){u c kR e dR )e ijk , (6) 

°f = (d c iL u jL )(u c kL e dL - d c kL v dL )e ijk , (7) 

Of = (d c iR u jR )(u c kR e dR )e ijk . (8) 

We can also list the operators relevant to strange final states ||: 

Of = (s C iR u iR)(v c kL e dL - d c kL is dL )e ijk , (9) 

Of = (s c iL u jL )(u c kR e dR )e ijk) (10) 

Of = (s c iL u jL )(u c kL e dL - d c kL u dL )e ijk , (11) 

Of = (s c iR u jR )(u c kR e dR )e ijk , (12) 

Of = (d c iR u jR ){s c kL u dL )e ijkl (13) 

Of = {d c lL u ]L ){s c kL v dL )e l0kl (14) 

where d denote the generation; e\ = e, e% = fi, v\ = v e and u 2 = v^. 

We are interested in the decay processes from the nucleon to one pseudoscalar meson: 
(proton,neutron)^(7r, K, rf) +(v, e + , /i + ). For these decay modes we can list the complete 
set of independent matrix elements in QCD employing the operators of eqs. (p|) — (|T1D : 

(ir°\e ijk (v? T CP R)L d j )P L u k \p), (15) 

(Tr + \e ijk (v? T CP R>L di)P L d k \p), (16) 

(K°\e ljk (u lT CP R , L si)P L u k \p), (17) 

(K + \e i3k (u* T CP R , L si)P L d k \p), (18) 

(K + \e ijk (v* T CP R>L (P)P L s k \p), (19) 

(K°\e ljk (u iT CP R , L si)P L d k \n), (20) 

( V \e ijk (u iT CP R , L d j )P L u k \p), (21) 

where we assume SU(2) isospin symmetry m u = m d and use the relations 

(PS\O lr \N) = (PS\O rl \N), (22) 

(PS\Orr\N) = (PS\O ll \N), (23) 
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due to the parity invariance. All we have to calculate in lattice QCD are these 14 matrix 
elements. Other matrix elements are obtained through the exchange of the up and down 
quarks, under which the nucleon and PS meson states transform as 



\p) -> -\n), \n) -> -\p), (24) 

(7r + |-(7r-|, (A -> -<tt°|, (7r-|-(ir + |, (25) 

(K + \ -> -> (if + |, (26) 

(77| - (77|, (27) 

where there is no decay mode with the K° or K~ final state. 



B. Form factors in nucleon decay matrix elements 

Under the requirement of Lorentz and parity invariance, the matrix elements between the 
nucleon(iV) and the pseudoscalar(PS) meson in eqs. (P~5|) — (PT|) can have two form factors: 

(PS(p)\Of\N^(k)) = P L (W (q 2 ) - W q {q 2 )i4) (28) 

where of represents the three-quark operator projected to the left-handed chiral state, 
denotes the Dirac spinor for nucleon with either the up (s = 1) or down (s = 2) spin state 
, and q 2 is the momentum squared of the out-going antilepton. The contribution of the W q 
term in eq. (|28|) is negligible in the physical decay amplitude, because its contribution is of 
the order of the lepton mass mi after the multiplication with antilepton spinor. However, 
since the relative magnitude of the two form factors Wq and W q is a priori not known, we 
have to disentangle these two form factors in the lattice QCD calculation. Hereafter we refer 
to Wq and W q as relevant and irrelevant form factor respectively. 

In the lattice calculation, k = is chosen for the nucleon spatial momentum and p = 
k — q^Oior the PS meson. In this case the Dirac structure of the right hand side in eq. ( |2"5| ) 
is given by 
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{w Q -w q i4)u 



Wo - iq^Wg -W q q ■ a 
W q q-a W + iq 4 W q 



u 



W + (m N - \Jm 2 PS + f)W q 

-W q p- a W 



W q p- a 



it 



(29) 



where Wq — W q i$ is expressed by a 2 x 2 block notation; a are the Pauli matrices, and 



u 



(1,0,0,0) or (0,1,0,0). It is important to observe that the upper components of 
(Wo — W q i£)vS 8 ' are linear combinations of the relevant and irrelevant form factors, while 
the lower components contain only the irrelevant one. Therefore, we can extract the relevant 
form factor Wq from the upper components by subtracting the contribution of the irrelevant 
form factor W q with the use of the lower components. 

The need for the separation of the contribution of the irrelevant form factor was not 
recognized in the previous studies with the direct method [§]||. The values found in these 
studies correspond to W — iq^W q instead of W - We examine how much this affects the 



estimate of the matrix elements in Sec. [VI. 

Let us add several technical comments: (i) The separation procedure described above 
cannot be applied to the case of p = k = because of vanishing lower components, (ii) 
Another possible choice of momenta for disentangling the relevant and irrelevant form factors 



is given by k ^ and p = 0. In this case, however, we cannot achieve — q 2 



mf. 



C. Calculational methods 

The nucleon decay matrix elements of eq. fl28|) are calculated with two methods, which we 
refer to as the direct and the indirect method. The former is to extract the matrix elements 
from the three-point function of the nucleon, the PS meson and the baryon number violating 
operator. The latter is to estimate them with the aid of chiral lagrangian, where we have 
two unknown parameters to be determined by the lattice QCD calculation. 

In the direct method we calculate the following ratio of the hadron three-point function 
divided by the two-point functions: 
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p E^e^^--)(J P5 (^,f)ag 7 (x,t)Jv, s (o)) z z 

e i P-( £ '- f )( J PS (f', t')J PS {x, t)) Y.x(Jn,s{Xi ^J'nA )) V ''' S ' V X 



1 (PS(p)|dfjjV< s >(fc = d)> t / »t>0. (30) 



L x L y L z x v 

Here C^ 7 denotes the renormalized operator in the naive dimensional regularization(NDR) 
with the MS subtraction scheme, and 7 and s are spinor indices; we can specify the spin 
state of the initial nucleon at rest by choosing s = 1 or 2. L x L y L z is the spatial volume of 
lattice in lattice units. The amplitudes Zps and are given by 



(PS(p)\Jp S (0,0)\0) = yZps, (31) 
(0|J w , s (0,0)|iV( s ')(0)) = Jz^u<f'\ (32) 



which can be obtained from the two-point functions 

£ e*l*-*>(Jp S (it, t')J PS (x, t)>, (33) 

^(Jjv iS (f, ^7^(0,0)). (34) 

We move the baryon number violating operator 0]° 7 in terms of t between the nucleon 

source placed at t — and the PS meson sink fixed at some t' well separated from t = 0. 

We list all the local interpolating fields for the PS meson and the nucleon required to 
calculate the independent matrix elements of eqs. Ql5|)— fl2"l~l): 

J T o (x, t) = -j= (u(x, t)j 5 u(x, t) - d(x, t)j 5 d(x, £)) , (35) 

J n +(x, t) = d(x, t)7 5 w(a;, t), (36) 

Jk°(x, t) = s(x, t)^^d{x^ t), (37) 

J K +(x,t) = s(x,t)-f 5 u(x,t), (38) 

Jr){x, t) = — j= (u(x, t)7 5 -u(x, t) + d(x, t)7 5 c?(x, t) — 2s(x, t)7 5 s(x, £)) , (39) 

J P;S {x, t) = e ljk (u iT (x, t)C-f 5 d j (x, tfj u k s (x, t), (40) 

J n ,s{x, t) = e ijk (u iT (x, t)C-f 5 d j (x, t)) d k s (x, t). (41) 
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We also prepare smeared operators for the nucleon source to overlap with the lowest energy 
state dominant ly: 

■W*) = E ^(x)^m(z)e ljk (u lT (x,t)C l5 di(y,t)) u k s (z,t), (42) 

x,y,z 

J'nM = E *(aWy)*(2)ey* (u lT (x,t)C l5 d^y,t)) d k s (z,t), (43) 



x,y,z 



where the measured quark wave function in the pion is employed for the smearing factor 
which is obtained by 



T,y{d{x, t)7 5 u(0, t)u(y, 0)7 5 d(y, 0)) 



i > (44) 



Eff(d(0, t)j 5 u(0, t)u(y, 0) l5 d(y, 0)> 
with configurations fixed to the Coulomb gauge. Although there is no reason to assume that 
the wave functions for the three quarks in the proton is well described by the quark wave 
function in the pion, the smeared sources of eqs. (|42]) and (|43|) work effectively (see Sec. |V|). 



In the renormalization of the baryon number violating operators on the lattice, the ex- 
plicit chiral symmetry breaking in the Wilson quark action causes mixing between operators 
with different chiralities. In eqs. (|T5|) — (f2~T|) we find two types of operators in terms of chi- 
ralities: 

Orl = e ijk {i,\ T CP R ^ 2 )P L <iPl (45) 
O ll = e l]k ^\ T CP L iPi)P L iPt (46) 

where ^1,2,3 represent the quark fields. Their mixing structures under perturbative renor- 
malization up to one-loop level are given by |5| 

O&iii) = Z(a s , pa)6%{a) + ^-Z mix O^(a) - ^Z' mix O^ L (a) , (47) 
Olt(v) = Z(a s , va)6%{a) + ^Z mix O^(a) + ^Z' mix O^ L (a), (48) 

where the overall factor Z(a s ,[ia) has the form 

Z(a s , H = 1 + ^ (41n(H + A 9 ) , (49) 

with fi the renormalization scale, and the additional operator 1 l is defined by 
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0^ L = e ijk {^\ CwMPl-y^I (50) 

Employing the MS subtraction scheme with the naive dimensional regularization for the 
continuum theory, we have reevaluated the finite constants and found 

A ¥ = -34.11 for NDR, (51) 
Z m ix — 3.21, (52) 

Z' mix = -0.803, (53) 

where the errors are ±1 in the last digit. The value of A^ depends on the renormalization 
scheme in the continuum, while Z mix and Z' mix are independent. We present the integral 
for A^ in the dimensional reduction(DRED) scheme and those for Z mix and Z' mix in Ap- 
pendix, where we give a detailed description for the one-loop perturbative calculation of the 



renormalization factors. With the use of the KLM normalization of quark fields [10] and the 



tadpole improvement [II], the overall renormalization factor of eq. (f49|) is rewritten as 

3 

Z{oc s , no) = (l- — ) 1 + ^ (iln(fxa) + A $ + -vr x 5.457^ . (54) 
Here K c is the critical hopping parameter at which the pion mass vanishes. We use 

1 1 - 5.457a s /4 (55) 



in Ref. |12| for the perturbative estimate of K c . 

Let us turn to the indirect method. The baryon number violating operators constructed 
in chiral lagrangian contains two unknown coefficients a and f3 defined by 



(0\e ijk (u^ CP R di)P L u k \p^) = olP l vS'\ (56) 
(0\e ijk (v* T CP L di)P L u k \p^) = f3P L u^ s \ (57) 



where operators are renormalized in the NDR scheme with the use of the renormalization 
factors of eqs. (|5"2"D— (f>§). These matrix elements are obtained from the two-point functions: 



Y J ^k{u tl CP R ^)P L u\x,t)J' p ^)) 

^ 3 (J P;S {x,t)J^{0)) 
(0\e ijk (u i CP RtL d j )P L u k \pW) t > 0. (58) 
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Incorporating the a and (3 values determined by the lattice calculation in the tree-level results 
of chiral lagrangian, we can evaluates the nucleon decay matrix elements of eqs. ( Dp — (pT|) . 



III. TREE-LEVEL RESULTS OF CHIRAL LAGRANGIAN 

For some nucleon decay matrix elements, tree-level results of chiral lagrangian have 
already been given in Refs. [|8],[13|, which are obtained with the use of the on-shell condition 
of the out-going leptons: — q 2 = mf and i^Vi = raiV\. In our lattice calculations, however, 
the lepton momentum is generally off the mass shell. Hence we need to understand the q 
dependence for an extrapolation of the matrix elements to the physical point. In this section 
we present the tree- level results for all the independent matrix elements in eqs. fllpp — (|21~D 
with the explicit expressions of q dependences. 

We first define the chiral lagrangian for baryon-meson strong interactions following the 

notation of Ref. ||. The PS meson and baryon fields are given by 

\ 



( J-ttO + -L n 

V2 71 + V6 71 

K~ 

( J_yo 1 J_ao 



71 

1 ^0 



73" ^Te 11 



7T U + 

K° 



K+ 
K° 



(59) 



B 



V 



Lyo 1 J_ ao 

wo 



V 
n 



(60) 



In terms of <fi we define the 3x3 special unitary matrices: 

'2i(j)\ 



exp 



(61) 
(62) 



where / is the pion decay constant. Under SU(3)iXSU(3)ij the meson and baryon fields 
transform as 



B -> UBU\ 
12 



(63) 
(64) 



where L is an element of SU(3)l and R is an element of SU(3)#; U is defined through the 
transformation properties of £ : 



The lowest order of the SU(3)lxSU(3)# invariant chiral lagrangian is given by 



(65) 



P 



Co = ±-Tr(d^)(d^) + TrB(j^ + M B )B 

o 



(66) 



on the Euclidean space-time. Quark mass contributions can be included by adding the 
symmetry-breaking term 



-v 3 Tr(^M q + M q E) 

- ai TrB(?M q ? + £M q £)B - a 2 TrBB(?M q ? + £M q £) 
-hTiB^M^ - ^M q i)B - b 2 TrB l5 B(^M q ^, -£M q £), 



(67) 



where 



M n 



\ 



m u 
m d 
m s j 
related to the meson mass by 



V 



(68) 



The parameter v is 

f 2 ml±fi _ _ Pm 2 K ± Pm 2 K0 

Mm,, + m s ) \(m u ■ ///,/ • 4m, 



4(m u + m d ) A(m u + m s ) 



3pml 



(69) 

Experimental results for the semileptonic baryon decays give F = 0.47 and D = 0.80 |I3 . 
The symmetry-breaking parameters a\ and a 2 are estimated from mass splittings among 
the octet baryons; b\ and b 2 , on the other hand, are not well determined since they do not 
contribute to the baryon masses. The parameters v, a± and a 2 have no contribution to the 
tree-level results for the nucleon decay matrix elements. 
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Let us consider the construction of the operators of eqs. (j5|) — (|i4|) , which are written 
in the quark fields, with the meson and baryon fields. The operators transform under 
SU(3) L xSU(3) i? , as 



(3,3) :Of,6f,6® 



d > 

(3,3): Of, 6f, 
(8,l):Of,af,(5f, 
(1,8) :Of,6f. 



(70) 
(71) 
(72) 
(73) 



These transformation properties are realized by £5£ G (3,3), G (3,3) G (8,1) 

and £^B£ G (1, 8), with which we can express the operators of eqs. (|5|)— (PH) as 



of 




v c dL TrF'£B L , 


(74) 


of 


= ae c dB Ti^B R e, 




(75) 


of 


= (3 (e c dL Ti^B L e - 


- VhTrPtBzg) , 


(76) 


of 


= i3e c dR Tr^B R ^ 




(77) 


of 


= a (e^Tr^B^ - 


- PhTtftBLt) , 


(78) 


of 


= ae c dR Ti^B R e, 




(79) 


of 


= P (e c dL Ti^B L e ■ 


- VhTrPtBzg) , 


(80) 


of 


= Pe c dR Tr^B R £, 




(81) 


of 


= ap c dL TrP'ZB L Z, 




(82) 


of 






(83) 



where a and /?, which are already defined in eqs. (|5q ) and (|57|) , are unknown coefficients 
associated with the (3, 3) and (3, 3) operators and the (8, 1) and (1, 8) operators respectively; 
T, J 7 ', J 7 , T 1 and T" are projection matrices in the flavor space, 



^ ^ 





1 



^ 0^ 




1 



(84) 



14 



/ 0^ 



\ 



-10 




^0 0^ 



0-10 




^0 ^ 




1 



(85) 



We can now apply the chiral lagrangian Cq + L\ and the baryon number violating op- 
erators of eqs. (|71|) — (|S3|) to calculating the nucleon decay matrix elements. Expanding the 
lagrangian and the operators in terms of the meson and baryon fields, we obtain the following 
tree-level results for the independent matrix elements of eqs. (p~5|) — (pl|): 

1 D + F—q 2 + m 2 N 4bi m u m N 



(ir \(ud R )u L \p) = aP L u p 



(n°\(ud L )u L \p) 



y\(ud R )d L \p) 



(lT + \(ud L )d L \p) 



(K°\(us R )u L \p) 



(K°\(us L )u L \p) 



y/2f y/2f ~q 2 - 



m 



N 



V2f ~q 2 - 



m 



N 



-aP L idu f 



D + F 2m N 
\/2f -q 2 - rr 



46i 



N 



y/2f ~q 2 



m 



N 



(3P L u p 
-pP L idu p 
aP L u p 
-aP L idu p 

-f3P L idu p 
aP L u p 
-aP L i<ju p 
(3P L u p 



D + F -q 2 + m 



V2f 



D 



V2f 
F 2m N 



q 2 - m 2 N 



46i m u m^ 
y/2f —q 2 — m 



N 



4bx 



m,, 



V2f 



f 

D + F 



m N 

-q 2 - m 2 N 



V2f ~q 2 



m 



N 



D + F-q 2 + m 2 N 2b x (m u + m d )m N 



2m 



N 



m 



N 



f —q 2 — m 2 
2bi m u + m d 
f ~q 2 
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where we use (PS\(ijjiijj2 R L )tp3 L \N) as a shortened form of (PS^^ 1 C P RL ^ 2 )P L ^\N)- 
q dependences of the matrix elements are retained without applying the on-shell condition of 
the out-going leptons — q 2 = mf, itfvi = m[V[. These expressions are considerably simplified 
if we employ the approximations of m s ~ m A = m B , m u ^ d <C m s , m s /m B <C m N /m B , 
h,2 ~ 0(1) and -q 2 < m NtB - 
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A. Data sets 



Our calculation is carried out with the Wilson quark action in quenched QCD at (3 = 6.0 
on a 28 2 x 48 x 80 lattice. Gauge configurations are generated with the single plaque- 
tte action separated by 2000 pseudo heat-bath sweeps. We employ 20 configurations for 
the measurement of the quark wave function in the pion, which is used for the nucleon 
smeared source, after the thermalization of 22000 sweeps, and then analyzed the next 100 
configurations for the calculation of the nucleon decay matrix elements. The four hopping 
parameters K = 0.15620, 0.15568, 0.15516 and 0.15464 are adopted such that the physical 
point for the K meson can be interpolated. The critical hopping parameter K c = 0.15714(1) 
is determined by extrapolating the results of m 2 at the four hopping parameters linearly 
in 1/2K to m\ = 0. The p meson mass at the chiral limit is used to determine the in- 
verse lattice spacing a -1 = 2.30(4)GeV with m p = 770MeV as input. The strange quark 
mass m s a = 0.0464(16) (K s = 0.15488(7)), which is estimated from the experimental ratio 
m K /m p = 0.644, is in the middle of K = 0.15516 and K = 0.15464. 



B. Calculational procedure 

Our calculations are carried out in three steps. We first measure the quark wave function 
in the pion for each hopping parameter using the ratio of eq. ((0). For this purpose we 



prepare gauge configurations fixed to the Coulomb gauge except the t = time slice. On 
these configurations the pion correlation functions in eq. ([14]) are constructed employing 
the quark propagators solved with wall sources at the t = time slice where the Dirichlet 
boundary condition is imposed in the time direction. We note that the non-local pion 
sources in the t = time slice cancel out in the average over gauge configurations. Figure [I] 
shows the results of ^(l^l) measured at t = 18 for the heaviest (K = 0.15464) and lightest 
(K = 0.15620) hopping parameters. We use the central values of ^/(|x|) for the smeared 



nucleon sources of eqs. (f42"D and 
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In the second step we calculate various two-point functions required to determine hadron 
masses, ^Zp$, y/Z^, a and (3. We extract the PS meson masses and the amplitudes y/Zps 
from the correlation functions of eq. fl33|) where we employ the set of quark propagators 
solved with the sources of e lp ' x at the t = 29 time slice without gauge fixing. The nucleon 
masses are determined from the smeared-local correlation function J2x(JN,s(x,t)J' N s (0)), 
fixing gauge configurations on the t = time slice to the Coulomb gauge. The amplitudes 
y/Z~N are evaluated by fitting the local-local correlation function of eq. fl34j) to an exponential 
form with the nucleon mass fixed. It is straightforward to calculate the a and (3 parameters 



with the use of the ratio of eq. (58 



Finally we calculate the ratio of eq. (p0|), where the baryon number violating operator is 
moved between the nucleon source and the PS meson sink. Gauge configurations on the t = 
time slice are fixed to the Coulomb gauge to employ the smeared source for the nucleon. For 
the calculation of the three-point function in the ratio, we use the source method to insert 
the pion fields at t = 29 into the quark propagators solved with the t — smeared source 



| T5| . We should note that calculation of the p — > rj matrix elements of eq. requires 
the disconnected diagrams in terms of the quark lines, which cannot be calculated by the 
source method. Although these diagrams could have contributions to the matrix elements 
in the non-degenerate case of the up, down and strange quark masses, we neglect them in 
this paper. Four spatial momenta pa = (0, 0, 0), (vr/14, 0, 0), (0, 7r/14, 0) and (0, 0, 7r/24) are 
imposed on the PS meson in the final state. For the p ^ cases we distinguish the strange 
quark mass from the up and down quark mass by providing different hopping parameters 
for rrii and m 2 in Fig. |2|. As explained in Sec. |II B| , we cannot disentangle the relevant form 
factor from the irrelevant one in the case of the PS meson at rest, where we take only the 
degenerate quark mass mi = m 2 . 

From the tree-level expressions of chiral lagrangian for the nucleon decay matrix elements 
in eqs. (|36D — (|99"D , we can assume that the form factors obtained from the ratio of eq. ( |30| ) 
are functions of q 2 , m x and m 2 , where the mi and m 2 dependences could appear through the 
baryon masses, the pion decay constant and the a, /3, F, D, &i j2 parameters, To interpolate 
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the form factors to the q 2 = point, where the charged lepton masses are negligible (see 
Sec. |VT|), we employ the following fitting function 

c + c x ■ {-q 2 ) + c 2 • {-q 2 f + c 3 -m 1 + c 4 - m 2 . (114) 
We extrapolate mi and m-i to the chiral limit for the matrix elements of eqs. (0), flT5| ) and 



u21|), while rri2 is interpolated to the physical strange quark mass with mi taken to the chiral 
limit for the matrix elements of eqs. fllTD — (|20|) . 

To calculate the perturbative renormalization factors, we determine the strong coupling 
constant at the scale 1/a and n/a in the MS scheme. We first define the coupling constant 
a P [[L(| from the expectation value of the plaquette P = (TrC/p}/3: 

47T 

- InP = — a P (3.40/a) [1 - 1.19a P ] . (115) 



The conversion from ap to the MS coupling constant is made by 

2 2 
a) 1 H — ap + 0.95a P 

7T 



<*Ms(3.40/a) = ap(es3.40/ 



(116) 



The values of a^(l/a) and a^g(7r/a) are obtained by two-loop renormalization group run- 
ning starting from aj^g(3.40/a). 

We estimate errors by the single elimination jackknife procedure for all measured quan- 
tities. 



V. RESULTS FOR a AND (3 PARAMETERS 

In this section we present the results for hadron masses, y/Zps, sfZ^ii a an d (3 which are 
obtained from the two-point functions. In Fig. |3] we plot effective masses of the PS meson 
for the case of \p\a = at K = 0.15464 and \p\a = 7r/14 at K = 0.15620; the statistical 
errors are best controlled in the former and worst in the latter. We observe plateaus beyond 
t ~ 13 for both cases. The horizontal lines denote the fitted values of the PS meson masses 
with an error of one standard deviation obtained by a global fit of the two-point function of 
eq. (|33D with the function 
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ZpS ( e -m PS (t-2%) + e -m PS (T-t+29)\ ^ (117) 

2mps ^ ' ' 

where the fitting range is chosen to be 13 < t < 22 after taking account of the time reversal 
symmetry t — 29 <-> T — t + 29. This fitting procedure also gives the amplitude Zp$. We 
tabulate the numerical values of mps in Tables p] and those of the PS meson energy for the 
case of | p\ 7^ in Table [TT]. 

Figure |] shows the nucleon effective masses obtained from the smeared-local correlation 
functions for the heaviest^ = 0.15464) and lightest^ = 0.15620) quark masses, which 
should be compared with Fig. [5] for the local-local correlation functions. We observe that the 
smeared source works effectively, dominantly overlapping with the lowest energy state. We 
extract the nucleon masses by fitting the smeared-local correlation functions to a exponential 
form with the fitting range 8 < t < 22. The fitted values are shown by the horizontal lines 
in Figs. (| and [| together with one standard deviation errors. The amplitudes Zjy defined in 
eq. (|32"D are obtained by a fit of the local-local correlation functions with the function 

Z N e~ mNt (118) 

over the range 13 < t < 22, where is fixed to be the value determined from the smeared- 
local correlation functions. We present the numerical values of for the four hopping 
parameters in Table |I[ 

The a and j3 parameters are extracted from a constant fit of the ratio of eq. (|5H), which 



is shown in Fig. ^] for the case of the lightest quark mass(i^ = 0.15620). The horizontal 
lines represent the fit with the fitting range chosen to be 8 < t < 22. The numerical 
values are given in Table |ITl[ Figure [7| illustrates quark mass dependences of the a and (3 



parameters. Applying linear fits to the data, we obtain a(NDR, 1/a) = — 0.015(l)GeV 3 and 
/5(NDR, 1/a) = 0.014(l)GeV 3 in the chiral limit with the use of a" 1 = 2.30(4)GeV. 

Let us compare our results for the a and (3 parameters with the previous estimates. We 
summarize the previous lattice results in Table [IV| together with the simulation parameters. 
In Refs. [||,U the lattice cut-off scale a -1 was determined by the nucleon mass. The nucleon 
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mass results employed |j!8,19| are, however, quite heavy compared to those of more recent 



high statistical calculations pO| , |2T| : m N a = 1.11(10) JT8| compared to m N a = 0.756(19) |2( 
in the chiral limit at (3 = 5.7, and m^a = 0.64(11) (19j compared to m^a = 0.461(9) |?T| in 
the chiral limit at (3 — 6.0. To avoid this large uncertainty, we employ a™ 1 determined by 
the p meson mass to obtain the a and (3 parameters in physical units in Table [TV|. 

In phenomeno logical GUT model analyses of the nucleon decays, the values |a| = \(3\ = 
0.003GeV 3 |2!| are conservatively taken as these are the smallest estimate among various 
QCD model calculations 0. A trend one observes in Table [IV| is that the previous lattice 
calculations indicated values of these parameters considerably larger than the minimum 
model estimate above. Our results, significantly improved over the previous ones due to 
the use of higher statistics, larger spatial size, lighter quark masses and smaller lattice 
spacing, have confirmed this trend: the values we obtained are about five times larger than 
101 = 0.003GeV 3 . 



VI. RESULTS FOR NUCLEON DECAY MATRIX ELEMENTS 

We now turn to the calculation of the nucleon matrix elements with the direct method. 
In Fig. | we show time dependences of R(t, t' = 29) with \p\a = ir/lA for the matrix element 
{K°\eijk{u lT CPRd?)PL,u k \p) in the case of the heaviest quark mass(i^ = 0.15464) and the 
lightest one(-?r = 0.15620). The results of constant fits are represented by the sets of three 
horizontal lines. We choose the fitting range to be 8 < t < 16 for all the matrix elements 
of eqs. (|15|) — (0) such that the excited state contaminations in the nucleon and PS meson 
states observed in Figs.[3] and £| can be avoided simultaneously. 

Figures [15] show —q 2 a 2 dependences of the relevant form factors Wo(q 2 ) in the indepen- 
dent nucleon decay matrix elements in eqs. (^) — (^TJ), where the operators are renormalized 
with the NDR scheme at fi = 1/a. The values of —q 2 a 2 are enumerated in Table [TT] as 
a function of the quark masses and the spatial momentum p. In Fig. ^| (a) we also 
plotted the combination W — iq^Wg for comparison, which is obtained by following the 
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method in Ref. |3J]. The magnitude of Wo(q 2 ) is more than two times larger than that of 
Wo(q 2 ) — ig^Wq {q 2 )- The relevant form factors at —q 2 a 2 = 0(open circles) in Figs. |9|— 
are obtained by fitting the data employing the function of eq. ( |114| ), where we find that 
the charged lepton masses m 2 a 2 = 4.9 x 10~ 8 and irv^a 2 = 2.1 x 10 -3 are negligible in the 
current numerical statistics. We plot the function c + c\ ■ (—q 2 ) + c 2 • (— q 2 ) 2 employing the 
fitting results of Co, c\ and c 2 in Figs. |9|, |10| and [15], and Co + c\ ■ (—q 2 ) + c 2 • (— q 2 ) 2 + c±- m s 



with the fitting results of Co, ci, c 2 and C4 in Figs Jill— |14|. We observe that the signs of Co 
and c\ are consistent with the predictions of chiral lagrangian in eqs. (|100|) — (|113|) for all 



the matrix elements, while the signs of c 2 show disagreement in some matrix elements. The 
coefficients c 2 , however, are poorly determined compared to c and c\. The fitting results 
for W (q 2 = 0) are presented in Table [V|. 

In Fig. [16] we compare the nucleon decay matrix elements obtained by the direct method 
with those by the indirect one using the tree-level results of chiral lagrangian (squares), 
where we employ the expressions of eqs. fllOOQ — ( |113| ) with a(NDR, l/o) = — 0.015(l)GeV 3 , 
/?(NDR,l/a) = 0.014(l)GeV 3 , / w = 0.131GeV, m N = 0.94GeV, m B = 1.15GeV, D = 0.80 
and F = 0.47. We observe that the two set of results are roughly comparable. This leads 
us to consider that the large discrepancy between the results of the two methods found in 
Refs. @,[| is mainly due to the neglect of the W q (q 2 ) term in eq. (|28f) . 

It is also intriguing to compare our results with the tree-level predictions of chiral la- 
grangian with |a| = \/3\ = 0.003GeV 3 (crosses) that is the smallest estimate among various 
QCD model calculations. Our results with the direct method are 3 — 5 times larger than 
the smallest estimates except (r]\(udji)uL\p) . Hence they are expected to give stronger con- 
straints on the parameters of GUT models. 

Finally, let us discuss the soft pion limit of the nucleon decay matrix elements. The tree- 
level result of the chiral lagrangian for (tt \{udR)UL\v) m ec L- ®) shows that the combination 
of form factors Wo(q 2 ) —iq4,W g (q 2 ) converges to a finite value of a/(\/2f) in the soft pion limit 
Pn — ¥ (— q 2 = m%), whereas each of Wo and W q diverges. In Fig. [T7| we plot Wq — iq^Wq 
for the matrix element (7r°\(udji)uL\p) as a function of —q 2 a 2 . In this case the results for the 

24 



pion at rest are also included. To extrapolate the data to the point —q 2 a 2 = m 2 N a 2 (dashed 
vertical line), we employ the fitting function 

Co + ci • (-g 2 ) + c 2 • {—iqi) + c 3 • mi + c 4 • m 2 . (119) 

Solid line denotes c + C\ • (—q 2 ) with the fitting results of c and C\. We also draw c + C\ • 
(— q 2 ) + C2 • (— iqi) + C3 -mi + C4 -m,2 (dotted lines) choosing the four cases of — ig4 = mjv — m n 
with mi = m 2 . The extrapolated value (open circle) at the point —q 2 a 2 = m 2 N a 2 is consistent 
with the result of a/(y/2f) (triangle). We observe similar situations for {^\{udL)ui\p) and 
{n + \{ud R)L )d L \p). 

VII. CONCLUSIONS 

In this article we have reported progress in the lattice study of the nucleon decay ma- 
trix elements. In order to enable a model-independent analysis of the nucleon decay, we 
have extracted the form factors of all the independent matrix elements relevant for the 
(proton,neutron)^(7r, K, rf) +(P, e + , fj, + ) decay processes without invoking chiral lagrangian. 

We have also pointed out the necessity of separating out the contribution of an irrelevant 
form factor in lattice calculations for a correct estimate of the matrix elements at the physical 
point. With this separation, the matrix elements obtained from the three-point functions 
are roughly comparable with the tree-level predictions of chiral lagrangian with the a and (5 
parameters determined on the same lattice. The magnitude of the matrix elements, however, 
are 3 to 5 times larger than those with the smallest estimate of a and (3 among various QCD 
model calculations. Our results would stimulate phenomenological interests as the larger 
values of the nucleon decay matrix elements can give more stringent constraints on GUT 
models. 

The ultimate goal of lattice QCD calculations of the nucleon decay matrix elements is to 
determine the matrix elements precisely with control over possible systematic errors. Major 
systematic errors conceivably affecting our present results are the scaling violations and the 
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quenching effects. The former can be investigated by repeating the simulation at several 
lattice spacings; the latter is eliminated once configurations are generated with dynamical 
quarks, where it is straightforward to apply our method. We leave these points to future 
studies. 



One of us (Y.K.) thanks C. Bernard for useful discussions. This work is supported by the 
Supercomputer Project No. 45 (FY1999) of High Energy Accelerator Research Organization 
(KEK), and also in part by the Grants-in-Aid of the Ministry of Education (Nos. 09304029, 
10640246, 10640248, 10740107, 10740125, 11640294, 11740162). K-I.I is supported by the 
JSPS Research Fellowship. 



The perturbative renormalization factors for the baryon number violating operators Orl 
and Oll, which are defined in eqs. ( pT/| ) and ([48]), have already been calculated in Ref. || 
employing the DRED scheme for the continuum theory. However, the authors of Ref. f{| 
present only the numerical results for Z, Z mix and Z' mix . We consider that it would be 
instructive to demonstrate the calculation of the renormalization factors in detail. 

We first rewrite the operators Orl and Oll as 



where tp c = ip T C is a charge conjugated field of if). The continuum and Wilson quark actions 
for the charge conjugated field ip c are obtained from those for ip with the replacement of 
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(120) 




(121) 



igT A 



^g(T A ) T ^ 



(122) 
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where T A (A = 1, ... ,8) are generators of color SU(3) group. This implies the modification 
of the Feynman rule of the quark-gluon vertex for the ip c field. 

We illustrate the relevant one-loop diagrams in Fig. [18]: (a) the quark self energy and 
(b)-(d) the three types of vertex corrections. We calculate these diagrams in the Feynman 
gauge for massless quarks with vanishing external momenta. The infrared divergences are 
regularized by introducing the fictitious small mass A in the gluon propagator: 

1 



n ab cont 



(123) 
(124) 



Ja ^4E a sin 2 (W2) + A2- 
We should note that the infrared behavior of the theory should be independent of the 
ultraviolet regularization schemes. The infrared divergent contributions in the one-loop 
diagrams, which emerges as the InA 2 terms, are supposed to cancel in the renormalization 
factors relating the continuum and lattice operators. 

Up to the one-loop level the inverse quark propagator and the vertex functions are written 
in the following form: 



G- 1 ( P ,A) = ^(l-^S«(A)), 



Arl,llW = Pr,l ®Pl + ll (A), 



(125) 
(126) 



where the superscript (i) refers to the i-th loop level. LL represents the sum of con- 
tributions from the three diagrams in Figs. |TE|(b) — (d). The continuum results for £W and 
A rl,ll are given by 



S«(A) 



— 7 + ln|47r 



In 



aS, ll (A) = 4P R)L ®P l 
in the NDR scheme and 



- 7 + ln|47r| + In 



A 2 



+ 



(127) 
(128) 



E«(A) 



7 + In 4?r + In 



A 2 



+ 



AgjU(A) = APrj. ® P L 



7 + ln|47r| ) + In 



+ 1 



(129) 
(130) 
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in the DRED scheme, where the reduced space-time dimension D is parameterized by e as 
D = 4 — e, e > 0. The pole term (2/e — 7 + ln|47r|) should be eliminated in the MS subtraction 
scheme. The corresponding lattice results for and A^| LL are 



S (1) (A) = -^ 



A<£ (A) = 41n 
2 



Ag(A) = 41n 
2 



4 


In 


A 2 






- 1 


3 




A 2 





+ ; 



2 ^ d 4 k 



l+r 2 



A 2 (4A! + A 2 ) 



A 4 + r 2 Ai(2- Ai) 



1 



1 +r z 



A 2 (4A X + A 2 ) 2 V 2 
-6(A 2 - k 2 " ~ 1 



A : A 4 - A 4 - A 5 + 



1 



2(4Ax + A 2 ) 



(fc 2 + A 2 ) 2 



(131) 



A=0 



A 2 



A 2 



d 4 k 



Pr®Pl 



(2vr) 4 

+2A 4 - 2A 5 + 6r 2 AiA 4 + 8r 4 A^ 
P L ®P L /A 



A|(4Ai + A 2 ) V 2 

4 A3 



^(4r 2 A 2 -A 4 ) + ^(4-(l-r 2 )A 1 ) 



^(4r 2 A 2 - A 4 ) - ^(4 - (1 - r 2 )A x ) + 2r 2 A x A 4 ) 



A 2 (4A X + A 2 ) 



-9(A 2 - k 2 ) 



6 



fc 2 (fc 2 + A 2 ) 



(132) 



A=0 



A 2 



+ 



A 2 

3 (4?r) Lite 



d A k 



A 2 (4Ai + A 2 ) 



(^(4r 2 A?-A 4 ) + ^(4-(l-r 2 )A 1 ) 



+2A 4 - 2A 5 + 6r 2 AiA 4 + 8r 4 A?) 



R 



(^(4r 2 A 2 - A 4 ) - ^(4 - (1 - r 2 )A x ) + 2r 2 A x A 4 ) 



A 2 (4A X + A 2 ) 



-9(A 2 - k 2 ) 



6 



A; 2 (A; 2 + A 2 ) 



(133) 



A=0 



where r denotes the Wilson parameter and Ai, A 2 , A 4 , A 5 and A 6 are given by 

Ax = ]T sin 2 (k a /2), (134) 
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A 4 = ^sin 2 (fc a ), 

a 

A 5 = ^sin 2 (fc a )sin 2 (£; Q /2), 

a 

A 2 = A 4 + 4r 2 A 2 , 
A 6 = (l + r 2 )Ax-4. 



(135) 
(136) 
(137) 
(138) 



The counter terms in proportion to 6(A 2 — k 2 ), which have the same infrared singularities as 
the lattice integrands, are introduced to pick out the analytical expressions of the infrared 
divergent contributions. The hyper-sphere radius A does not exceed 7T. With the use of 
eqs. ( |127| ) — ( |130| ) and ( |131|) — ( |"133|) , we obtain the expression for the renormalization constant 
As in eq. 



NDR 



A 



DRED 



+ 



~-41n|A| 



-(4tt 



d 4 k 



- w (2tt) 4 
1 /1 + r 2 



1 + r* 



8A 2 A 



2 \ 2 
1 /A 6 



6AiA 2 



2AiA 2 V 8 
AiA 4 - A 4 - A 5 
A 2 



A 4 + r 2 Ai(2 - A x ) 
1 

+ 



(4r 2 A 2 



A 4 ) + =^(4 



4A a 



+2A 4 - 2A 5 + 6r 2 AiA 4 + 8r 4 A? 
2 



-^(a 2 - e 



Zmix = (47r) 2 

zl x = w 2 r 

J —IT 



(139) 



The mixing coefficients Z mix and Z' mix in eqs.(|47D and ( (4"8| ) are expressed by 

- (SeA^f (t* 4 ^ 2 " A <> " T< 4 " " r2 > A '' + 2r2A ' A ' ) ' < 140 > 



- d*k 1 (r 2 



-7T (2vr) 4 6A| V 4 



— A 4 - 4r^Ax + r l A\ . 



(141) 



We evaluate numerical values of A« DR , Z mix and Z' mix with r = 1 using the Monte Carlo 



integration routine BASES 



which are already presented in Sec. |11 Q Our results for 



Z mix and Z' mix are consistent with those in Ref. [f|, while we observe a slight deviation 



beyond the statistical error of the numerical integration for A^ 



DRED 
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FIG. 1. Quark wave function in the pion normalized by the value at the origin for (a) 
K = 0.15464 and (b) K = 0.15620. \x\ is distance between two quarks. 
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FIG. 2. Quark flow diagram for the nucleon decay three point function with the mass assign- 
ment. Filled circles denote the local operators and shaded rectangular is for the smeared operator. 
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FIG. 3. (a) Effective mass for the pion with \p\a = at K = 0.15464 and (b) effective energy 
for the pion with \p\a = 7r/14 at K = 0.15620. The pion correlation functions consist of the local 
sink and the wall source without gauge fixing. Solid lines denote the fitting results with an error 
band of one standard deviation obtained by global fits of the pion propagators. 
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FIG. 4. Effective mass for the nucleon with the smeared source at (a) K = 0.15464 and (b) 
K = 0.15620. Solid lines denote the fitting results with an error band of one standard deviation 
obtained by global fits of the nucleon smeared-local propagators. 
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FIG. 5. Effective mass for the nucleon with the local source at (a) K = 0.15464 and (b) 
K = 0.15620. Solid lines denote the fitting results with an error band of one standard deviation 
obtained by global fits of the nucleon smeared-local propagators. 
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FIG. 6. Ratio R a/3 (t) for (a) a and (b) (3 parameters at K = 0.15620. Solid lines denote the 
fitting results with an error band of one standard deviation. 
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FIG. 7. Chiral extrapolations of (a) a and (b) (3 parameters. Solid lines denote linear fits. 
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FIG. 8. Ratio R(t,tf = 29) for the relevant form factor in {n Q \{ud R )u L \p) at (a) K = 0.15464 
and (b) K = 0.15620. Solid lines denote the fitting results with an error band of one standard 
deviation. 
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FIG. 12. Same as Fig. y for (a) (K + \(us R )d L \p) and (b) (K + \{us L )d L \ 
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FIG. 14. Same as Fig. O for (a) (K°\(us R )d L \n) and (b) (K°\(us L )d L \n) 
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FIG. 18. One-loop diagrams for (a) quark self energy and (b)— (d) vertex corrections for the 
three-quark operator, p denotes a external quark momentum and i, j and k at the ends of quark 
lines label color indices. 
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TABLES 



TABLE I. Hadron masses at (3 = 6.0 in quenched QCD. 



K 



0.15464 
0.15516 
0.15568 
0.15620 
0.157136(12) 



0.3209(8) 
0.2843(9) 
0.2436(9) 
0.1957(11) 



0.4350(22) 
0.4135(26) 
0.3925(33) 
0.3723(45) 
0.3346(56) 



0.6674(38) 
0.6253(43) 
0.5814(52) 
0.5359(69) 
0.4607(89) 
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TABLE II. Four-momentum transfers from the nucleon at rest to the pseudoscalar meson. Eps 



is the energy 


of the pse 


:udoscalar meson 


with spatial momentum p. 










\p\a = 


\p\a - 


7T/ Z4 


\p\a 


= tt/14 


K for m\ 


K for 777,2 


-q 2 a 2 


tips • a 


— a 


E PS ■ a 


-q 2 a 2 


0.15464 


0.15464 


0.1200(28) 


U.o4o / (1U J 


U.UoO ( (zo J 


0.3912(12) 


0.0259(22) 




0.15516 




n Qono/ 1 1 ^ 
U.ooUz(ll J 


u.uyoo^z ( J 


0.3767(13) 


0.0341(24) 




0.15568 




U.oiol(lz J 


U.lUo4(zy J 


0.3619(15) 


0.0429(25) 




0.15620 




U.zy0o(14J 


n 1019/Q1\ 

U.lzlo(ol J 


0.3469(21) 


0.0523(28) 


0.15516 


0.15464 




U.ooU4(lU J 


U.Ubyo(z/ J 


0.3769(13) 


0.0113(23) 




0.15516 


0.1163(31) 


n Q1 Qoh 1 ^ 
U.oloz(liJ 


U.U8Uo(zy J 


0.3620(14) 


0.0190(25) 




0.15568 




U.zyOo(lo J 


U.(jyio(oU J 


0.3469(15) 


0.0272(26) 




0.15620 




U.Z (00(10 J 


U.lU40(oZ J 


0.3316(22) 


0.0359(28) 


0.15568 


0.15464 




U.oioo(li J 


U.Uo40(zy J 


0.3626(14) 


-0.0025(24) 




0.15516 




u.zyo / (iz j 


U.Ub40(ol J 


0.3473(16) 


0.0045(26) 




0.15568 


0.1141(36) 


U.zYoo(14J 


U.U ID I (66 ) 


0.3317(19) 


0.0120(28) 




0.15620 




0.2567(16) 


0.0883(35) 


0.3159(25) 


0.0202(31) 


0.15620 


0.15464 




0.2969(13) 


0.0400(33) 


0.3480(18) 


-0.0150(27) 




0.15516 




0.2777(14) 


0.0495(36) 


0.3322(19) 


-0.0089(29) 




0.15568 




0.2575(15) 


0.0604(39) 


0.3163(23) 


-0.0021(32) 




0.15620 


0.1158(47) 


0.2356(18) 


0.0730(42) 


0.3003(31) 


0.0052(37) 
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TABLE III. Results for a and parameters as a function of quark mass. Operators are 
renormalized at the scale \i in the NDR scheme. 



era 3 0a 3 



K 


fi = 1/a 


fi = n/a 


[i = 1/a 


fi = ir/a 


0.15464 


-0.00180(6) 


-0.00207(7) 


0.00179(6) 


0.00205(7) 


0.15516 


-0.00166(7) 


-0.00191(8) 


0.00165(7) 


0.00189(8) 


0.15568 


-0.00155(7) 


-0.00177(8) 


0.00152(7) 


0.00174(8) 


0.15620 


-0.00148(9) 


-0.00170(11) 


0.00143(9) 


0.00164(10) 


0.157136(12) 


-0.00125(11) 


-0.00144(13) 


0.00119(11) 


0.00137(12) 



TABLE IV. Comparison of a and parameters in lattice QCD. All calculations are done with 



the Wilson quark action in the quenched approximation. 


Lattice cutoff a 


1 is determined from 


m p . Quark mass is 


defined by m q = 


(1/2K - 1/2K C ) ■ a' 1 . 






Ref. 


Hara et al. g[F 


^ Bowler et al. @,|l8| 


Gavela et al. [^ 


19] This work 


Lattice size 


16 3 x 48 


8 3 x 16 


10 2 x 20 x 40 28 2 x 48 x 80 


No. config. 


15 


32 


30 


100 


a- 1 [GeV] 


1.81(6) 


1.45(9) 


2.2(4) 


2.30(4) 


Spatial size [fm 3 ] 


(1.7) 3 


(1.1)" 


(0.9) 2 x 1.8 


(2.4) 2 x 4.1 


Quark mass [MeV] 


109^^696 


184^m^477 


82^m g ^223 


44^m 9 ^118 


a [GeV 3 ] 




\a\ ~ 0.065 


\a\= 0.019(2) a = -0.015(1) 


[GeV 3 ] 


\/3\ = 0.029(6) 


\/3\ ~ 0.050 




= 0.014(1) 


Renorm. scheme 




Pauli-Villars 


DRED 


NDR 


Renorm. scale 




n = 85 GeV 


fi = 1/a 


fi = 1/a 
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TABLE V. Results for relevant form factors in the independent nucleon decay matrix elements 
of eqs. (|i~5|) — (^l"l) . Operators are renormalized at the scale ji in the NDR scheme. 





fj/ 1 / Li 




jJL — 7T / d 






laince uniLo 


uc V 


laince LiiiiLs 


vje v 


(tt \(udji)uL,\p) 


— U.UZDo^ol J 




n n98Q(''?fi\ 
— u.uzoy [oo j 


n i ^fi q"\ 
— u. ioo^ iy ) 


\7T \{UaL)UL\P) 


U.UZ^fcZl OZ 1 


D 1 98f 1 7"! 


n n978('^7 N \ 

U.UZ / O^O / J 


U. 14e ( ^ZU j 


/ 'Tr 1 1 { 1 1 /-/ -r* \ /-/ r 1 nr\\ 

\7r \{ua,R)aL\P) 


— U.UOO 1 \Qd J 


— u. ioy ^z^t j 


— u.u^tuy^oi j 


n 91 (\( r >7\ 


[IT \{UaL)UL\P) 




nisi ( r >A\ 


u.uoy^t^oz j 


fl 908^98^ 


[ft \{USR)UL\p) 


u.uiyz^zu ) 


n 1 n9fi i 


n n9i i( r > r >\ 

U.UZlo^ZZ J 


U.llO^lZ J 


I Tf®\ (n, a i Xn r ln\ 

[ft \{USl)Ul\P) 


n nnsoi' 1 9^ 
u.uuoy^iz ) 




u.uuyz^io ) 


u.u4o i ^oy j 


[ft \{USR)CIL\P) 


— U.UlUU^lOj 


— u.uozy ^t>y ) 


n ni 1 s^i ^\ 


— U.UDZ4H i yj 


[ft \{USl)CLL\P) 


n nnssi' 1 1 


U.U4eDD^0o j 


n m ri^fi ^ 

U.UlUO^lO J 


u.uooo^oy j 


[ft \\uaR)SL\P) 


— U.UZ4rO^Z / J 


— U. lol ^14t ) 


D 0989/^1 "\ 




(K + \(ud L )s L \p) 

\ IV iv / J-J \F 1 


0.0268(28) 


0.142(15) 


0.0304(32) 


0.161(17) 


(K°\(us R )d L \n) 


0.0090(12) 


0.0476(63) 


0.0094(13) 


0.0497(69) 


{K°\(us L )d L \n) 


0.0179(18) 


0.0947(95) 


0.0198(20) 


0.105(11) 


{r]\(ud R )u L \p} 


0.0007(11) 


0.0037(58) 


0.0004(12) 


0.0021(63) 


(rj\(ud L )u L \p) 


0.0191(25) 


0.101(13) 


0.0212(28) 


0.112(15) 
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